Possibility of narrow resonances in nucleon-nucleon channels 
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Compound states manifest themselves as bound states, resonances, or primitives, and their char- 
acter is determined by their interaction with the continuum. If the interaction experiences a pertur- 
bation, a compound state can change its manifestation. Phase analysis of nucleon-nucleon scattering 
indicates the existence of primitives in the ^^i, ^^o, and ^Po channels. Electromagnetic interaction 
can shift primitives from the unitary cut, turning them into narrow resonances. We evaluate this 
effect on the ^So proton-proton scattering channel in the framework of the Simonov-Dyson model. 
We show that electromagnetic interaction turns a primitve with a mass of 2000 MeV into a dibaryon 
resonance of approximately the same mass and a width of 260 keV. Narrow resonances of a similar 
nature may occur in other nucleon-nucleon channels. Experimental confirmation of the existence of 
narrow resonances would have important implications for the theory of nucleon-nucleon interaction. 



PACS numbers: ll.80.Et, 12.39.Ba, 13.40.Dk, 13.75.Cs 



Analytical properties of scattering amplitudes were ex- 
tensively studied in the 1950s, and for some time, the 
conditions of analyticity and unitarity were considered 
likely to be sufficient for the full recovery of scattering 
amplitudes. Low []] found the dispersion relation for 
scattering amplitudes, which takes into account analytic- 
ity and unitarity. Solving this equation, Castillejo, Dalitz 
and Dyson observed an ambiguity, which is now called 
the CDD poles. The physical meaning of these poles was 
clarified by Dyson . Using a version of the Lee model 
he showed that the CDD poles correspond to bound 
states and resonances (for a review, see Q). 

S'-matrices with a finite number of poles, which arc 
mainly used in the phenomenological parameterizations 
of scattering amplitudes, correspond to zero-range singu- 
lar potentials. The character of the interaction in such 
systems is determined by analyzing the behavior of the 
scattering phase rather than the potential: if the phase 
increases with energy, the interaction is an attraction; if 
it decreases, the effect of interaction is a bound state or 
repulsion. According to the Breit-Wigner formula, iso- 
lated resonances drive the phase shift up. In the absence 
of bound states, the systems discussed in Refs. [H-Q cor- 
respond to phases that increase with energy. 

In nucleon-nucleon interaction, the only bound state, 
the deuteron, arises in the ^Si channel, whereas all other 
two-nucleon channels are free of bound states. 

Conversely, the nucleon-nucleon phase shifts decrease 
with increasing energy and provide evidence for repul- 
sion. The Dyson model 3 was extended by Simonov @ 
for the description of internucleon forces. In the general- 
ized model, both attraction and repulsion may dominate, 
and the phase shifts behave bidirectionally. The repul- 
sion has been included by weakening the condition of 
strict positivity of the imaginary part of the D function, 
thereby allowing zeros of the D function to appear on the 
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unitary cut 7]. 



Numerous experimental searches for exotic multiquark 
states did not give decisive results. In the late '70s Jaffe 
and Low Q proposed an experimental method to iden- 
tify exotic hadrons using a special formalism, called the 
P-matrix. The application of the P-matrix formalism 
to meson-meson scattering revealed P-matrix poles that 
roughly correspond to the four-quark states predicted 
earlier by quark models P . States that show up as poles 
of the P matrix, rather than the S matrix, are called 
" primitives" . They correspond to zeros of the D function 
on the unitary cut, modify the Low scattering equation, 
and generate the CDD poles [S]. 

The Simonov-Dyson model provides a dynamical 
framework of the P-matrix formalism. It has been suc- 
cessfully applied to describe the nucleon-nucleon inter- 
action dominated by repulsion at short distances jol-lisj. 
In the Simonov-Dyson model, primitives are the objects 
that produce repulsion. 

The state of the art that prevailed in the early '80s 
can be characterized as follows: The Dyson model de- 
scribes systems dominated by attraction where bound 
states and resonances may exist. The Simonov-Dyson 
model describes systems with both attraction and repul- 
sion where, in addition to bound states and resonances, 
primitives come into play. 

The conventional approach to the description of 
nucleon-nucleon interaction is based on the Yukawa 
meson-exchange mechanism. As a result of the devel- 
opments of the '70s and early '80s, a second approach 
has emerged, which ascribes nucleon-nucleon interaction 
to the s-channel exchange of 6(7-primitives (for a review, 
see [13). 

The mechanism of t-channel meson exchange is studied 
in detail and is widely used to analyze particular effects. 
The s-channel exchange mechanism 15] has solid the- 
oretical background and is also largely self-sufficient in 
modeling nucleon-nucleon interaction. 
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The above two mechanisms may appear mutuaUy ex- 
clusive. To decide which of them is more appropriate, it 
is necessary to find a prediction that could draw a clear 
distinction between the alternatives. The task is not sim- 
ple, because they both successfully describe a wide range 
of experimental data. 

However, these mechanisms can also be dual to each 
other in the sense in which the s-channel resonance ex- 
change is dual to the i-channel resonance exchange in the 
Veneziano model. In our case, this would mean that the 
i-channcl meson exchange parameterizes the s-channel 
exchange of 6g-primitives, and vice versa. Simultaneous 
treatment of the t- and s-exchanges would represent in 
this case double counting. 

There is, moreover, a third possibility. The hybrid Lee 
model can be reduced to the Simonov-Dyson model [9|, 
so the latter can be regarded as a phenomenological real- 
ization of the former. In the hybrid Lee model, the inter- 
nuclcon forces are generated by the s-channel exchange 
of 6(7-compound states. 

We wish to clarify which of the aforementioned mech- 
anisms is realized in nature, and whether there is a du- 
ality between the t-channel exchange of mesons and the 
s-channel exchange of primitives (6g-compound states). 

The existence of primitives imposes constraints on the 
interaction parameters. A characteristic feature of the 
Simonov-Dyson model is that small perturbations of the 
interaction with the continuum shift primitives from the 
unitary cut and, in the general case, turn them into res- 
onances. 

Although primitives do not exist in the OBE models as 
elementary objects, the behavior of the phase shift asso- 
ciated with the repulsion permits a natural interpretation 
in terms of primitives. Any intersection of the phase shift 
of level (5(s) = mod (tt) with a negative slope can be 
interpreted as a zero of the D function and in turn as a 
primitive @ . The status of primitives in the OBE models 
is roughly the same as the status of resonances and bound 
states. Small perturbations of the t-channel exchange pa- 
rameters move primitives along the unitary cut. If there 
is a duality, the zeros of the D function do not leave the 
unitary cut, so primitives remain primitives. The duality 
condition thereby imposes constraints on the permissible 
perturbations of the parameters of the Simonov-Dyson 
model, which determine, in particular, compound state 
masses and their coupling with the continuum. 

A similar situation occurs in the hybrid Lee model. 
The Simonov-Dyson model ensures the existence of sta- 
ble primitives, if it arose as a reduction of the hybrid 
Lee model. In this case, arbitrary perturbations of 
the parameters of the hybrid Lee model are allowed, 
whereas the corresponding variations of parameters of 
the Simonov-Dyson model are those that keep primitives 
on the unitary cut. 

Some discrimination between the aforementioned al- 
ternatives can be made, if one considers scattering in an 
external field. In the OBE models and in the case of 
duality, the scattering cross section will not be qualita- 



tively changed. Conversely, if there is no duality and the 
s-channel exchange mechanism operates as assumed in 
the Simonov-Dyson model we can expect the appearance 
of a resonance in the place where we earlier observed a 
primitive. Because the mass of the resonance is approx- 
imately known (it corresponds to zero phase shift), the 
problem reduces to the estimation of the width of the 
resonance. In the hybrid model of Lee, primitives are 
rigidly connected with the unitary cut, so the resonance 
will not occur. 

In fact, we can proceed without external fields. 
A meaningful example of the perturbation of strong- 
interaction parameters of the model is given by the elec- 
tromagnetic interaction. In this case, the condition for 
the existence of primitives is violated as a result of 
electromagnetic mass shifts of the nucleons and the (Sq- 
compound states, as well as modification of the imaginary 
part of the D function, which additionally receives the 
Gamow-Sommerfeld factor for the proton-proton chan- 
nel. Inclusion of the electromagnetic effects allows the 
determination of the width of the resonance, as well as 
the (small) shift of the resonance mass relative to the 
primitive mass. 

We therefore consider the possibility that the displace- 
ment of primitives is associated with electromagnetic in- 
teraction. This takes place in all nucleon-nucleon chan- 
nels because of the electromagnetic mass shifts and other 
electromagnetic corrections, including the initial- and 
final-state Coulomb interaction. For numerical estimates, 
we use the Simonov-Dyson model [l6j . 

The D function can be written as 



where 
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The indices (5 label 6g-compound states with masses M^, 
the coupling constants parametrize the interaction 
with the continuum, and / describes the four-fermion 
contact interaction. $2(3) — nkj^fs is the phase space 
volume of two nucleons, where k is the nucleon momen- 
tum in the center-of-mass frame. J^(s) is the form factor 
arising from the vertices of compound states and nucle- 
ons and the four-fermion interaction vertices: 



T{s) = - 
so 



1/4 



sin(fc6) 
kb 



Co(s). 
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Here, sq is threshold value of s and b is the efi^ective 
interaction radius. The Gamow-Sommerfeld factor [13, 
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exp(27r?7) — 1 
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where rj = a/j/k and fj, = m/2 is the reduced proton 
mass, accounts for the initial- and final-state Coulomb 
interaction of the 5'-wave protons. In the pn and nn 



channels, Co(fc) = 1. 

In the proton-proton channel, the dispersion integral 
of Eq. ([3|) can be evaluated to give 



J 



2k% 
sin(fc6) 
kb ' 



(2Si(2fc&) sin(2fc&) -I- (Ci(2fcfo) -I- Ci(-2fc6)) cos(2fc6) - 2C - \n{~4:k^b^)) 



ikb 



sinh(b/n) ^„b/„ 
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where k = 2mb/'K and C = 0.5772 ... is the Euler con- 
stant. Here units are chosen in which Cc = l/(a^) — 1, 
later we will restore the dimension. In the expression 
(|n]), a branch of the logarithm ln(l) — enters; in the 
complex fc-plane, the cut extends from —oo to 0, so 
ln(— 1 ± iO) — iiTT. In the physical upper half of the 
complex fc-plane, the self-energy operator n(s) is an an- 
alytical function. On the real s-axis below sq or, equiva- 
lently, on the imaginary half-axis SRfc = 0, 5fc > 0, n(s) 
is a real function. 

A^^(s) is assumed to have a single pole at s = M^, 
corresponding to a Gg-compound state. The contribu- 
tion to A~^(s) of the contact four-fermion interaction is 
described by a constant term. As a result. 
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So 



M-^' 7' 
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where the coupling constant entering 

Eq. ([2]) and rp = Sn'^/b'^ is the residue of the free P 
matrix. The value of Cp measures the strength of inter- 
action of the compound state with the continuum. The 
residue of k~^A~^(s) is parameterized as follows 



So 
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In the pn and nn channels, the D function has the correct 
analytical properties, provided ^ £ (0,1). The phase 
shifts are well reproduced for ^ = 0.9. Here we restricted 
ourselves to verification of the fact that, in the channel 
at ^ = 0.9, the analytical properties of the D function 
are correct. 

The value of 7 is determined by the scattering length. 
In the pp channel, we use the series expansion around 
fc = for 3?fc > fii: 



kCl{s) cot 6{k) 



2h{k) = - + -k^ 
a 2 



(8) 



This expansion allows one to draw a link between the 
low-energy parameters of the model and the phenomeno- 
logical constants that characterize the low-energy proton- 
proton scattering. The following equation serves as the 
definition of an analytical function h{k): 

2^{1 + y) + ik + ln(-fc2) = 2h{k) + ikC^is), (9) 



with ip{z) being the digamma function. 

The strong-interaction phase shift can be found from 
the strong-interaction 5* matrix modified by the initial- 
and final-state Coulomb interaction 



S'{k) 



D{s-iO) 
D{s + iO)' 



(10) 



We use the empirical value of the scattering length a 
to fix 7: 
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m^2 
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6'"C(m), (11) 



where C(''^) is the Riemann zeta function. Once the 
model parameters are determined, we arrive at a pre- 
diction for the effective range 



3 ^ '^'^ 6 {A'P ~ 4to2) 



7, 



(12) 
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, 2"'+2(m-l)(m + 6) 1 
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In the ^5*0 pn channel with Qc = 00, one has r = 2.1 fm 
0, which should be compared with r — 2.8 fm [22] in 
the OBE models. 

In bare strong interaction, the real and imaginary parts 
of the D function simultaneously vanish at some point on 
the real fc-axis. The zeros of 3_D(s), on the other hand, 
are the zeros of J-{s), the lowest-lying one is located at 

M 



n/bo. 



(13) 



In the parameterization ([7]), diD{s) also vanishes at this 
point, whereas the compound state mass is equal to the 
total energy of protons in the center-of-mass frame: 



Mo = 2Jnybl 



(14) 



In strong plus electromagnetic interaction (strong in- 
teraction + QED), the zero of the D function is shifted 
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to the complex plane. The effect is related to the change 
of the model parameters Bq ^ b — bo + A6, mg m — 
Too + Ato, etc., and to the Coulomb interaction of the 
protons. Let us discuss the change of parameters. 

Electromagnetic mass splitting of hadrons is well stud- 
ied. QED corrections to hadron masses are usually at- 
tributed to the Coulomb interaction between the quarks 
and their spin-spin interaction: 



AA^ = c < ^Cjej > 



-h < eiejCTiiTj > . (15) 



Here, it is assumed that all the quarks are in S'-wave 
states. The averaging is made over the color-spin-isospin 
wave function of the quarks. A reasonable description of 
the mass splitting in octet-baryon isomultiplets, account- 
ing for the mass difference of nonstrange quarks, can be 
obtained for c = 3.06 MeV and h = -1.35 MeV 

The average values of the spin-isospin operators and 
the corresponding electromagnetic mass shifts of the pro- 
ton, the neutron, and the 6g-compound states with quan- 
tum numbers of the ^Sq nucleon-nucleon channels are 
listed in Table ID The calculation of average values of 
the operators was carried out using fully antisymmetrized 
color-spin-isospin quark wave functions. The resonance 
is denoted by d* . 

For the nonstrange quarks, the MIT bag model gives 
b ^ R and M ^ B? , where R is the radius of the 6- 
quark bag, and thus the change of the parameter b for the 
inclusion of electromagnetic interaction is simply given 
by 



A& 



AM 
"3M ' 



(16) 



The values of the ^Sq scattering lengths, derived from 
the experimental data, are as follows [22|: 7.83 ±0.01 fni, 
23.75 ± 0.01 fm, and 16.4 ± 0.09 fm in the I3 = 1, 0, -1 
channels respectively. In our approach, the scattering 
length or, equivalently, the value of 7 is adjustable pa- 
rameter. We obtain a = 6.5 fm in the pp channel by 
fitting a maximum of the phase at T « 10 MeV. 

To compare the physical phase with that of the bare 
strong interaction, we must have an estimate of uq. We 
give some heuristic arguments in favor of the change in 
the scattering length being proportional to the electro- 
magnetic mass shift of the 6(7-compound state with quan- 
tum numbers of the channel. 

In the Born approximation, the scattering length is 
proportional to the averaged interaction potential. A 
part of this potential associated with the electromagnetic 
interaction leads to a difference of scattering lengths in 
different isospin channels. On the other hand, the elec- 
tromagnetic mass splitting of hadrons is also determined 
by averaging the electromagnetic potential, although of 
quarks rather than nucleons. In any case, the equation 
a ^ ao + C2AM, with oq = 14.45 fm and C2 = -1.45 
fm/MeV, gives scattering lengths a = 7.3 fm, 23.4 fm, 
and 16.6 fm in the I3 = 1,0,-1 channels, respectively. 



TABLE I: Electromagnetic mass shifts of tlie proton, the 
neutron, and the 6g-compound states d* with isospin projec- 
tions 1 3 = ±1,0. The third and fourth columns show average 
values of the spin-isospin operators. Electromagnetic mass 
shifts of the hadrons are shown in the last column (in MeV). 
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which agrees well with the empirical data. Thus, we ac- 
cept 



Aa = C2AM. 



(17) 



A slightly better fit of the proton-proton phase shift is 
obtained for a = 6.5 fm. Because a/uc ~ 1/3 (in dimen- 
sional units flc — 58 fm) , the electromagnetic corrections 
to low-energy parameters of the model are expected to 
be only moderately small. In the scattering lengths, the 
electromagnetic corrections ~ 1. 

Parameter ^ is also an adjustable parameter. We will 
accept A^ = 0. 

We thus have estimates of the proton mass toq in the 
bare strong interaction, as well as the electromagnetic 
mass corrections for nucleons and d*. The relative value 
of the correction to b is known also. 

The proton kinetic energy Ts = 244 MeV in the lab- 
oratory frame corresponds to a vanishing phase shift. 
In the center-of-mass frame, one finds the momentum 
of the protons kg = \J mTs /2 and their total energy 
Ms = \/2^^i7F±4^ = 1995 MeV. The phase shift van- 
ishes provided that the imaginary part of the D function 
vanishes; that is, it is necessary to have kgb — tt, so 
the interaction radius b is fixed. The total mass of the 
compound state M w Mg is known with an accuracy of 
0(AM), so we find from Eq. ((T6)) the interaction radius 
60 with electromagnetic interaction switched off. 

The condition (IT5|) holds for bare strong interac- 
tion. Using Eq. (1141) one may find a mass of the com- 
pound state with electromagnetic interaction switched 
off. Then, for a 7^ 0, the mass M = Mo + AM. 

In Fig. [1] the proton-proton scattering phase shift is 
plotted in strong interaction ± QED and in the bare 
strong interaction case. In Fig. [2J the S'-wave proton- 
proton cross section is plotted. The narrow resonance 
peak is associated with the complex root of the equation 



(18) 



D(M2-iM,r*) = 0. 



If the electromagnetic interaction is switched off, Eq. ([T? 
has a primitive root corresponding to F* = 0. When the 
electromagnetic interaction is switched on, a resonance 
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of a width 7^ appears in the neighborhood of the 
primitive. 

The fuh set of the model parameters is shown in Table 
im In the first row, the model parameters to which we 
previously attributed the subscript are given. 

The value of Ms is defined to within a few MeV. If we 
take Ms = 2000 MeV, the resonance will have a mass of 
= 2005.5 MeV and a width of = 250 keV. Because 
the position of the peak is not very well defined, in an 
experimental search for the resonance, one needs to scan 
a region around the zero value of the scattering phase 
with a resolution in -^i better than 100 keV. 

Phase analysis of nucleon-nucleon scattering did not 
reveal primitives in the "^Pi and ^Pi channels. In the 
■^Po channel, the data indicate the existence of a prim- 
itive with a mass of about 1970 MeV [l^. Arguments 
similar to those used for the S'-wave scattering allow one 
to conclude that electromagnetic interaction in the ^Pq 
proton-proton channel transforms the primitive to a res- 
onance with a mass of about 1970 MeV and a width of a 
few hundred keV. 




0.2 
T [GeV] 

FIG. 1: (color online) The ^So proton-proton scattering phase 
shift in radians modulo vr versus the proton kinetic energy in 
the laboratory frame. The smooth dashed curve 1 is the 
model prediction with electromagnetic interaction switched 
off. The solid broken curve 2 arises upon dressing the model 
parameters by electromagnetic interaction, i.e., in strong in- 
teraction -|- QED. The squares show the experimental data 



In nuclear matter, nucleons are in an effective mean 
field. Here we should expect significant renormalization 
of the parameters responsible for the dynamics. In the 
Simonov-Dyson model, we also should expect the trans- 
formation of the Og-primitives to the resonances, which 
significantly modifies the properties of nuclear matter. 
Dibaryon resonances are bosons. In cold nuclear mat- 
ter and with increasing chemical potential of nucleons, 
a Bose condensation will begin [2J]. This modifies the 
equation of state (EoS) of nuclear matter and modifies 
the properties of neutron stars under equilibrium p5| . 
Bose condensation of dibaryon resonances may also be 



the trigger for a phase transition to quark matter. The 
above scenario and its connection with other effects in 
nuclear physics and astrophysics are discussed in Ref. 

M. 
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FIG. 2: (color online) The ^5*0 proton-proton partial- wave 
cross section versus the proton kinetic energy in the labora- 
tory frame. Notations are the same as in Fig. [T] 



Astrophysical observations of neutron stars provide 
important constraints on the nuclear matter proper- 
ties above saturation density and, indirectly, on the in- 
medium modifications of hadrons. The problem of the 
possible conversion of primitives to resonances, however, 
can most reliably be settled in the laboratory. 

In the phase behavior of the nucleon-nucleon scatter- 
ing, there are clear features that can be attributed to 
primitives. We therefore believe that in bare strong in- 
teraction some compound states remain on the unitary 
cut and manifest themselves as primitives, rather than 
resonances or bound states. When external conditions, 
such as non-zero density of matter or finite temperature, 
become different or new interactions come into play, com- 
pound states can leave the unitary cut, turning into res- 
onances. These resonances can be observed in the labo- 
ratory as the usual Breit-Wigner peaks. 

There is no conclusive evidence for the existence of ex- 
otic 6q resonances. However, in the phase analysis of NN 
scattering 6q primitives are determined reliably. In the 
P-matrix formalism primitives are identified with exotic 
multiquark states. If primitives are shifted from the uni- 
tary cut under influence of perturbations, they can be 
observed as narrow resonances. 

In the OBE models and the hybrid Lee model primi- 
tives remain on the unitary cut. In the Simonov-Dyson 
model, in general, the primitives, if exist, are shifted from 
the unitary cut. Some preferences for these models can 
be done after the search results of the narrow dibaryons 
will become available. 

For an experimental search for narrow dibaryon reso- 
nances, one should have a beam of protons with a kinetic 
energy of T ~ 250 MeV and an energy spread below 



6 



TABLE II: Parameters of the Simonov-Dyson model for 
AIs ~ 1995 MeV with electromagnetic interaction switched 
off and on. a is the scattering length, r is the effective range, 
and b is the interaction radius (in fm) . m is the proton mass, 
M is the compound state mass, and M« and F* are the prim- 
itive and resonance masses and widths (in MeV). 
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100 keV. The energy resolution of the detector is not im- 
portant, so one can use scintillation detectors. In the 
CELSIUS accelerator at Uppsala, the beam momentum 
spread was a few times 10~^ before electron cooling and 
a few times 10~^ after electron cooling. Under such con- 
ditions, it is possible to measure ^/s with an accuracy of 
10 keV or better. The narrow width of d*++ (2000) is not 
an obstacle to its experimental search. 

Resonances of the same nature should exist in the pn 
and nn channels. Experiments with neutron beams are, 
however, more complicated. Also, e.g., deuterium cannot 
be used as a fixed target, because the narrow resonances 
are smeared out by the Fermi motion. We thus estimated 
the effect in the proton-proton channel, which is of inter- 
est for experiments with protons bombarding a hydrogen 
target. 



Suppose that, in proton-proton scattering, experimen- 
talists scanned the energy region up to a few MeV around 
the value where the phase shift is zero and the resonance 
was not found. In this case, we would have to conclude 
that the primitives are stable under perturbations and 
are strictly adhered to the unitary cut. 

Such stability is inherent in the OBE models and 
hybrid Lee model. Consequently, the Simonov-Dyson 
model would either be dual to the OBE models or be 
a phenomenological realization of the hybrid Lee model. 

Assume instead that, in place of the primitive, experi- 
mentalists found the d*"'""'"(2000) resonance with a width 
of about 260 keV. This would mean that the primitives 
are not strictly adhered to the unitary cut. This property 
holds only in the Simonov-Dyson model. 

Observation of the resonance will have the following 
consequences: i) mutual transformation of primitives 
and resonances will be confirmed, ii) universality of the 
Yukawa meson-exchange mechanism will be questioned 
on the experimental basis, iii) the region of validity of the 
hybrid Lee model will be reduced, and iv) the absence of 
duality between the i-channel meson exchange and the 
s-channel exchange of 6g-primitives will be proved. 

The author is grateful to M. G. Schepkin and Yu. A. Si- 
monov for useful discussions. This work is supported in 
part by RFBR grant No. 09-02-91341 and grant of Sci- 
entific Schools of Russian Federation No. 4568.2008.2. 
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